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Abstract 

Analytic expressions for the amplitude-dependent tune shift driven by the quadrupole mag- 
^ . net fringe field have been obtained. The magnitude of the effect is compared with other sources 

of non-linearity such as chromatic sextupoles, octupole errors of the main quadrupole mag- 

■ netic field and kinematic terms. A numerical example was the calculation on the lattice of the 
VEPP-4M e+e- collider. 

S' 

Ph! 1 Introduction 
o 

^ The main sources of non-linear detuning for storage rings and synchrotrons are the low-order mul- 
^ tipole components of the magnetic field (sextupole and octupole ones). Many authors investigated 
, the influence of these components on beam dynamics both through numerical simulation and ana- 
^1 lytically. The expressions estimating the non-linear betatron tune shift for such perturbations are 
r-| ■ well-known. 

, There are, however, other sources of the amplitude-tune shift that are able to influence signifi- 

cantly the beam behavior. For instance, the non-linearity of the fringe field was considered in [Ij for 
the low-energy ring (LER) of the PEP-II facility. It was shown that the main contribution (larger 
6o . than that of the regular chromatic sextupoles) was made by the fringe field of the quadrupole mag- 
^ ' nets. The kinematic effects and the fringe field of the bending magnets have no significant influence. 

■ For the i?-factory of KEK [2], it is pointed out that the main limits of the dynamical aperture are 
. determined by the kinematic terms of the interaction point drift space and by the fringe field of the 

O , quadrupoles adjacent to the interaction point. 

Q \ The studies are usually performed by computer simulation when the map over the quadrupole 

edge is written in terms of the Lee operators |3j or with the help of a second-order matrix formalism 

Below is made an attempt to obtain simple analytic expressions for the contribution of the 
quadrupole magnet fringe field to the non-linear tune shift. The formulae for the kinematic effects 
and ordinary octupole non-linear detuning are also presented for comparison. The octupole non- 
linear component is considered here as a small error of an ideal field inside the quadrupole lens. The 
contribution of perturbations of various kinds is compared on the basis of the VEPP-^M electron- 
positron storage ring. 

2 Hamiltonian 

We will consider a quadrupole magnet with an edge field drop (see, for instance, [5]). We will assume 
that an ideal quadrupole field inside the magnets is disturbed only by an octupole field error. Hence, 
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The transversal components of the quadrupole magnetic field are expressed as 



-- 

Bx = 

where Bi = (dB/dx) 



Bi{s)x + Ib3{s){x^ - 3xz^) - ^B'l{s){x^ + 3xz^) + 0(5), 
o 12 



B,{s)z + lB;{s)i3xh - z^) - ^B'{{s){z' + 3xh) + 0(5) 
12 



1 = 2=0 



is the field gradient, -B3 = {d^B/dx^ 



x=z=0 



(2.1) 

is the octupole component 



and B'{ = d^B/ds^ describes the contribution of the fringe field. Since the octupole field and edge 
field non-linearity are presented by the same power expansion series (but with different signs), the 
latter is referred to as pseudooctupole non-linearity sometimes. 

The Hamiltonian of the transversal motion of an on-energy relativistic electron in a magnetic 
field (EH) has the form of 0: 



B = \{pi +P2) + \ki{s){x'^ — z^) - the undisturbed linear part, 
+ \{p'i + pIY ^ kinematic part, 



-j^k'l{s){x^ + 6x^z^ - z^) 

-\k[{s)x'^zp^ 

+ ^A;3(s)(x4-6a;V + z^) 



the influence of the fringe field, 
the octupole component, 



(2.2) 



where ki{s) = Bi{s)/Bp and ^3(5) = B^{s)/Bp. It is considered here that fci > for the horizontally- 
focusing quadrupole magnet. 

The main aim of this work is to investigate the non-linearity of the fringe field. However, for 
the purpose of comparison and estimation, expressions for the kinematic terms and for the regular 
octupole component of the magnetic field will be also obtained. 

In order to find the amplitude-tune dependence in the first order of approximation, let us write 
down expressions (12. 2p in the "action-angle" variables [7] with the help of the generating function 



F{x,z,(t)xAx,s) 



tancp, — 



X 



2/?.(s) 



tan — 



(2.3) 



which assigns the following relations for the old and new variables: 

-^;|^(sin + ax{s) cos<^. 



X 

z 



2Jx(3x{s) cos( 



X ; Px 

2 J^(3-,{s) COS (f)^, Pz = - 



2Jz ( 

Ms) V 



sm( 



+ 0:2(5) cos I 



(2.4) 



where a and (3 are the Twiss parameters. 

Averaging the new Hamiltonian over all phase variables, (if(s))^ 
tion tune of the system according to 
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d 



ds. 



2tt I dJx 

we obtain the following first-order expressions for the amplitude-tune dependence: 

^^x CxxJx ~l~ CxzJz} 



and computing the oscilla- 

(2.5) 



(2.6) 



where the coefficients include the contribution of the kinematic effects, fringe field and regular oc- 
tupole perturbation: C* = O'^ -|- O^ -|- C°. Coefficients of each type have the following form: 
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Ci = ^/T,(s)T.(s)<is, (2.7) 

where 7 = (1 + a^)/ (3. 

b) Fringe field coefficients. 

= / K{s)f3l{s)ds, CI, = +^ / k';{s)(3!{s)ds, (2.8) 

= -^jp^{s){k'[{s)^{s)-Ak[{s)a,{s))ds. 

c) Regular octupole component coefficients. 

= ^jk,{s)Pl{s)ds,Cl, = ±-jk,{s)Pl{s)ds, (2.9) 



C:, = ^jk,{s)(3,{s)f3,{s)ds. 



For the sake of clearness and comparison, let us simplify general expressions (I2.7l) - fl2.9p basing 
on the evaluation of each type of perturbation in a cyclic accelerator. 

3 Kinematic effects 

The main contribution to integral expressions (12. 7p is made by the drift space with an extremely 
small value of the betatron function. The interaction region is best suited to consideration from such 
point of view (see, for instance, |2j). For a straight section with the length L, where behavior of the 
beta-functions is mirror-symmetric relative to the center, 7 = l//5o = const {(3q is the value at the 
center of the section), expressions fl2.7p take a simple form, convenient for estimation: 



4 Octupole errors 

If the magnet lattice of the accelerator is known and, thus, the behavior of the beta-functions along 
the beam trajectory is calculated, one can easily compute coefficients (12.90 numerically. However, 
we obtain the following simple estimation to be able to compare the influence of the octupole error 
of the central part of a quadrupole lens L long with the influence of the fringe fields of this lens. 

Let us assume, for the coefficient C°^, that Px is approximately constant inside the lens (fci ^ 1) 
and the octupole component k^ is constant over the length, then 

^ -r^k.PlL. (4.1) 

The quality g of a magnetic field is usually defined as the relative difference between the real 
and ideal values of the field at the radius a (we will regard this radius as the inscribed radius of the 



3 



lUiiowy iivdi 




(4.2) 



Substituting (14.21) into (14. ip . we obtain the estimation we will need later: 




(4.3) 



5 Fringe fields 

The non-linear betatron tune shift caused by the fringe fields of lenses is the principal subject of this 
work. In order to simplify expression (12. Sp . we need to introduce a model describing the behavior of 
the fringe field of the quadrupole lens. 

We have considered the following variants of model distribution of the quadrupole lens fringe 
field: the one based on M.Bassettis's formulae [Hi t9|, the model of "four lines" by G. Lee-Whiting 
[To] and the description of the magnetic field drop with two matched parabolas [H]. The first two 
models describe rather precisely the fringe field and its derivatives with respect to the longitudinal 
coordinate. 

For instance, for the "four-line" model, the distribution of the gradient and its first two derivatives 
(required for computation of (12.81) ) at the edge of the lens has the following form: 



where a is the inscribed radius of the aperture of the lens and Biq is the gradient on the lens axis in 
the central section. The approximation of the fringe field with the help of (15. ip is shown in Fig lS.H 
15.21 Unfortunately, it is rather complicated to use (15. ip since the field distribution is assumed infinite 
on each side of the axis (the edge corresponds to s = 0). So, a correct integration of (12.80 with due 
regard to the real behavior of the betatron function becomes difficult. 
At the same time, a simple model of two matched parabolas pLj 



yields a rather good approximation of the resulting expressions without cumbersome evaluation of 
numeric coefficients for the model (15. ip . The difference between the value of coefficients for models 
(15. ip and (15.21) will be shown below. 

FigEm 15.21 present the measured gradient drop of the quadrupole lens [T2J and the model ex- 
pression according to (15. ip and (15. 2p . The main notations used below are also given. The inscribed 
radius of the aperture of the lens a = 2.8 cm; the maximal gradient at the center of the lens Bio — 20 
T/m. It is worth noting that because of the quadrupole symmetry the end-pole chamfer, usually 
used for adjustment of the efficient length of the lens, has practically no influence on the fringe fleld 
pseudooctupole component deflned by B[{s). 

The fringe effects are concentrated in the vicinity of the longitudinal edge of the lens with a typical 
region length close to the lens inscribed radius: ~ ±a. We assume that the betatron functions vary 




(5.1) 




(5.2) 
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Experiment 
Model 1 
Model 2 




Figure 5.1: Fringe gradient drop: measured and computed, according to fl5.ll) (model 1) and according 
to dOD (model 2). 



Quadrupole fringe field B/'(s) 



Experiment 
Model 1 
Model 2 




Figure 5.2: Second derivative of the gradient with respect to the longitudinal coordinate. Notations 
are in Fig J5.1[ 
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PAs)-P..+PL{s-Se), (5.3) 

where and are the magnitudes of the beta-function and its derivative at the edge and e = 1, 2 
for the first and second edges. 

Let us substitute (15.31) and corresponding derivatives of the fringe field (15. 2p into (12.81) and 
perform the integration. Then we can obtain the following expressions describing the contribution 
of the fringe field of the quadrupole lens to the non-linear shift of the betatron tune: 

CI. = ^ho{^,f3'x,-f3xA.-/3.A+PxA2)^ 

CL = ^hoiP.A.-P.A^)^ (5.4) 

where kio is the focusing coefficient of the lens at the center. It is taken with its proper signs for 
the focusing and defocusing lenses. It is worth noting that the expressions obtained include only 
the parameters of the linear optics and do not include the lens aperture a. The complete value of 
the non-linear dependence of the tune on the amplitude is obtained by summing fl5.4p over all the 
quadrupole lenses of the accelerator. 

Usage of more accurate distribution of the fringe field (15.11) and integration of (12.81) close to the 
edge of the lens in the interval ±na gives the following correction in (15. 4p : kiQ should be replaced 
with 

"10 



2 (^2 + 1)3/2 ^2 + 1^' 

i.e. with n = 2 the difference is 7% and with n = 3 it is as small as 2%. 

Formulae (15.41) can be easily used numerically if the distribution of the betatron functions along 
the ring is available. However, to determine the general features of the influence of the quadrupole 
fringe field upon the non-linear shift of the betatron tune, it would be useful to simplify (15. 4p still 
further assuming that ^/ki^L < 1 for the quadrupole lens under consideration. Then (15. 4p can be 
expanded in series to the second order inclusive. 

^xx ~ ■J" ^10/^^0-^(1 ~ 7^10-^^), 

OTT 4 



1 .0 .0 1 



^ —kiAL{l + -k,oL'), (5.5) 

where /?o is the extreme (maximum or minimum) value of the corresponding betatron function 
achieved at the center of the lens. The last assumption is not essential and is made only because 
of the especially simple form of formulae (15. 5p in that case. In the case of arbitrary behavior of the 
betatron function (still having extremes inside the lens), in the brackets should be replaced with 
11 — I1I2 + I21 where li^2 is the distance from the lens edges to the extreme of the betatron function. 

One can see from (15. 5p that the frequency- amplitude dependence is always positive for the fringe 
field of the lens (at least, for the first-order term of the expansion), be this lens a focusing or a 
defocusing one. The last fact is illustrated well by Table iQl obtained numerically for LER pj. 

Let us compare the infiuence of the octupole error of the field inside the lens with that of the 
pseudo-octupole fringe component. From (14. 3 p and (15. 5p one can obtain: 
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Figure 6.1: Behavior of the lattice functions in the region of the VEPP-4M final focus. 

Substituting reasonable parameters of the quadrupole lens in (15.51) . one can see that the infiuence 
of the fringe field is comparable with and can even exceed the field errors in the inner area. 

6 Comparison of different contributions to the non-linear 
tune shift 

The magnet lattice of the electron-positron storage ring of VEPP-^M [13] was used as a model; the 
following comparisons were made for this machine: 

• The magnitudes of the non-linear shift of the betatron tune for perturbation of different kinds, 
including the fringe field, octupole error of field of the lenses, kinematic effects and sextupole 
magnets, compensating the natural chromaticity; 

• The tune dependence on the amplitude for the fringe field of the quadrupole lenses (obtained 
by the numerical modeling and in the analytic form, with the use of (15. 4p and (15.51) ). 

The octupole component of the field of the quadrupole lenses was taken from the results of the 
magnetic measurements |12j|. Its most significant value was in the defocusing lens of the final focus 
{EL2 in FigjnH). 

From the magnetic mapping, = 0.2 G/cm^ for this lens with an inscribed aperture diameter 
of 170 mm at an energy of 1.85 GeV, for which the calculations were performed. For the rest lenses 
of the ring this value is from 7 to 10 times less. 

Table 16.11 presents the comparative data for different sources of non-linearity for VEPP-^M 
It is seen that for the case of VEPP-^M the influence of the fringe fleld is relatively small. 
However, the corresponding values obtained for LER pQ by the numerical simulation exceed the 
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Kinematic term 




Figure 6.2: Cxx due to the kinematic effects for 1/2 of the VEPP-^M ring. The interaction point 
corresponds to s = 0. 



Table 6.1: Values of coefficients (EH) for VEPP-4M. 





C.., (m-i) 


a., (m-i) 


C.,, (m"i) 


Kinematic effects 


50 


7 


96 


Fringe field 


170 


300 


390 


Sextupole lenses 


88 


-1680 


-1660 


Octupole errors 


-510 


-340 


-280 
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Figure 6.3: C^x due to the fringe effects for 1/2 of the VEPP-^M ring. 



Table 6.2: Numerical simulation and analytical evaluation of the fringe effects 





C.,, (m-i) 




C.z, (m-i) 


Simulation 


182 


308 


403 


Evaluation bv 05.41) 


170 


300 


390 


Evaluation bv (15.51) 


170 


260 


350 



VEPP-4M values several times and their contribution to the non-linear shift is the most important 
among all kinds of perturbation. For instance, C^x from the sextupoles equals 800 m~^ while that 
from the fringe field effects is 1628 m~^. 

In case of kinematic effects, as expected, the main contribution is made by the interaction drift 
space (Fig J6.2l) . The quadrupole lenses ELI and EL2 of the final focus respond mainly to the fringe 



field non-linearity because the beta-functions here are from 110 to 120 m (Fig j6.3p . 

In order to compare the evaluation of the fringe field with the numerical simulation, the latter 
was carried out for the case of a "piecewise" approximation of the quadrupole lens field. In so doing, 
according to P, HI] , both the particle coordinate and the momentum change stepwise at the edge of 
the lens: 



Ax = —hoix^ + ^xz"^) 



^Px 



12 



2xzp^ -Px{x^ + z^] 



(6.1) 



For the vertical plane the substitutions x ^ z and fcio ^ — ^lo are made. The results of the simulation 
and analytical evaluation are compared in Table 16. 2[ 

It is seen that the analytic evaluation differs from the simulation by several percent. This can 
be explained by the accuracy of extraction of the tune shift at the Fourier discrete analysis of 1024 
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of the quadrupole lenses agrees rather satisfactorily with the results of the numerical simulation. 
Thus, one can obtain the result quickly from the linear magnetic structure of the cyclic accelerator, 
especially as few codes of the numerical tracking allow taking into account the influence of fringe 
non-linearity for the simulation. 
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